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Abstract 

We investigate zero sound in D-dimensional effective holographic theories, whose 
action is given by Einstein-Maxwell-Dilaton terms. The bulk spacetimes include 
both zero temperature backgrounds with anisotropic scaling symmetry and their 
near-extremal counterparts obtained in 1006.2124 [hep-th], while the massless charge 
carriers are described by probe D-branes. We discuss thermodynamics of the probe 
D-branes analytically. In particular, we clarify the conditions under which the spe- 
cific heat is linear in the temperature, which is a characteristic feature of Fermi 
liquids. We also compute the retarded Green's functions in the limit of low fre- 
quency and low momentum and find quasi-particle excitations in certain regime 
of the parameters. The retarded Green's functions are plotted at specific values 
of parameters in D = 4, where the specific heat is linear in the temperature and 
the quasi-particle excitation exists. We also calculate the AC conductivity in D- 
dimensions by-product. 
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1 Introduction 



The AdS/CFT correspondence [HE] has revealed the deep relations between gauge theo- 
ries and string theories and has provided powerful tools for understanding the dynamics 
of strongly coupled field theories in the dual gravity side. In recent years, this paradigm 
has been applied to investigate the properties of certain condensed matter systems [3]. 
The correspondence between gravity theories and condensed matter physics(sometimes is 
also named as AdS/CMT correspondence) has shed light on studying physics in the real 
world in the context of holography. 
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It is well known that in realistic condensed matter systems, the presence of a finite density 
of charge carriers is of great importance. According to the AdS/CFT correspondence, the 
dual bulk gravitational background should be charged black holes in asymptotically AdS 
spacetimes. The simplest example of such charged AdS black holes is Reissner-Nordstrom- 
AdS(RN-AdS) black hole, which has proven to be an efficient laboratory for studying 
the AdS/CMT correspondence. For instance, investigations of the fermionic two-point 
functions in this background indicated the existence of fermionic quasi-particles with 
non- Fermi liquid behavior [H El IE], while the AdS2 symmetry of the extremal RN-AdS 
black hole is crucial to the emergent scaling symmetry at zero temperature [TJ. Moreover, 
adding a charged scalar in such background leads to superconductivity [8j [9j [10] . 

A further step towards a holographic model-building of strongly-coupled systems at finite 
charge density is to consider the leading relevant (scalar) operator in the field theory side, 
whose bulk gravity theory is an Einstein-Maxwell-Dilaton system with a scalar potential. 
Such theories at zero charge density were analyzed in detail in recent years as they mimic 
certain essential properties of QCD [HI [121 [131 EH EH] • Solutions at finite charge density 
have been considered in [161 E3 EH1 ESJ [201 El] in the context of AdS/CMT correspondence. 

Recently a general framework for the discussion of the holographic dynamics of Einstein- 
Maxwell-Dilaton systems with a scalar potential was proposed in [22], which was a phe- 
nomenological approach based on the concept of Effective Holographic Theory (EHT). 
The minimal set of bulk fields contains the metric g^, the gauge field and the scalar 
(j) (dual to the relevant operator). <p appears in two scalar functions that enter the effec- 
tive action: the scalar potential and the non-minimal Maxwell coupling. They studied 
thermodynamics of certain exact solutions and computed the DC and AC conductivity. 
The main advantage of this EHT approach is that it permits a parametrization of large 
classes of IR dynamics and allows investigations on important observarables. However, it 
is not clear whether concrete EHTs can be embedded into string theories. For subsequent 
generalizations see [23] EH E5] E6J EH [28] ESJ [30]. 

On the other hand, strongly coupled quantum liquids play an important role in condensed 
matter physics, where quantum liquids mean translationally invariant systems at zero (or 
low) temperature and at finite density. By now there are two successful phenomenological 
theories of quantum liquids: Landau's Fermi-liquid theory and the theory of quantum 
Bose liquids, describing two different behaviors of a quantum liquid at low momenta and 
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temperatures. In particular, the specific heat of a Bose liquid at low temperature is 
proportional to T q in q spatial dimensions, while the specific heat of a Fermi liquid scales 
as T at low T, irrespective of the spatial dimensions. 

One may wonder if the newly developed techniques in AdS / CFT correspondence can help 
us understand the behavior of quantum liquids. In [31] the authors considered a class of 
gauge theories with fundamental fields whose holographic dual in the appropriate limit 
was given in terms of the Dirac-Born-Infeld (DBI) action in AdS space. They found 
that the specific heat~ T 2p in p spatial dimensions at low temperature and the system 
supported a sound mode at zero temperature, which was called "zero-temperature sound" . 
One interesting feature was that the "holographic zero sound" mode was almost identical 
to the zero sound in Fermi liquids: the real part of the dispersion relation was linear in 
momentum (u = qv) and the imaginary part had the same q 2 dependence predicted by 
Landau. The crucial difference was that the zero-temperature sound velocity coincided 
with the first-sound velocity, while generically the two velocities are not equal for a Fermi 
liquid. Such analysis was performed in the case of massive charge carriers in [32] and in 
the case of Sakai-Sugimoto model in [33]. The specific heat of general Dp/Dq systems 
was calculated in [31] and the specific heat of Lifshitz black holes was discussed in [35] 
and [36], while the zero sound was also investigated in [36]. 

In this paper we will study the low-temperature specific heat and the holographic zero 
sound in effective holographic theories. Here the bulk effective theory is D-dimensional 
Einstein gravity coupled to a Maxwell term with non-minimal coupling and a scalar. It 
was found in [37] that the theory admitted both extremal and near-extremal solutions 
with anisotropic scaling symmetry. We consider dynamics of probe D-branes in the above 
mentioned backgrounds and find that by appropriately fixing the parameters in the effec- 
tive theory, the specific heat can be proportional to T, resembling a Fermi liquid. We also 
compute the current-current retarded Green functions at low frequency and low momen- 
tum, and clarify the conditions when a quasi-particle excitation exists. Moreover, we also 
explore the possibility of observing the existence of Fermi surfaces in such a system by 
numerical methods. We find that although the system possesses some features of Fermi 
liquids, such as linear specific heat and zero sound excitation, we do not observe any 
characteristic structure in the wide range of k. In addition, the AC conductivity is also 
obtained by-product. 
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The rest of the paper is organized as follows: the exact solutions of the effective bulk 
theory will be reviewed in section 2 and the thermodynamics of massless charge carriers 
will be discussed in section 3. We shall calculate the correlation functions in section 4 and 
identify the quasi-particle behavior, while the existence of Fermi surfaces will be explored 
in section 5 via numerics. We will calculate the AC conductivity in section 6, including 
the zero density limit. Finally we will give a summary and discuss future directions. 

2 The solution 

In this section we will review the solutions obtained in [37], which can be seen as gen- 
eralizations of the four-dimensional near-extremal scaling solution discussed in [22]. In 
the beginning we consider the following action in .D- dimensions, without any reference to 
string theory or M/theory origin nor specifying the forms of the gauge coupling f(4>) and 
the scalar potential V(0) explicitly, 

S = -TT^yr- ! d D x^[R + f(<j ) )^^ + hdc^f + V(0)]. (2.1) 
lb7rG D J 2 

The resulting solutions are charged dilaton black holes, which have been investigated in 

the literature for a long period [381 EHJ EQJ HI]. Let us focus on solutions carrying electric 

charge only. The general configuration with planar symmetry can be written as follows 

dr 2 ^ 

ds * = _ u{r)d e + ^_ + v{r) Y^dxi 

^ ' 1=1 

(j> = 0(r), A t = A t (r), A r = Ai = 0. (2.2) 

After plugging in the scaling ansatz 

U(r) ~ r p , V(r) ~ r 7 (2.3) 

into the equations of motion, we can arrive at several constraints on the parameters and 
the scalar functions: 

• We require that /3 > 1, so that the extremal solutions have smooth connections to 
the finite temperature solutions; 

• The field equations indicate that /3 < 2 and < 7 < 2 and (3 = 2 when the scale 
invariance is restored. 
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• The equations of motion determine that the scalar field must take the form 

0(r) = C 2 logr + o , (2.4) 
and both /(</>) and V(4>) are constrained to be exponential in 0, power law in r. 

• Once we have fixed > 1, the metric must have /3 > 7, where saturation occurs for 
vanishing flux, e.g. in AdSo with /3 = 7 = 2. 

Subsequently, according to the constraints discussed above, we take the following forms 
of f(<f>) and V ((/>). 

f(4>)=e a4 >, V(4>) = -V e^. (2.5) 



Then we will consider charged dilaton black holes with a Liouville potential [12]. Now 
the scaling ansatz turns out to be 



„2 D-2 



0(r) = C 2 logr + o , A\r)= ^ ^ . (2.6) 

Since 0o an d C3 can be eliminated by rescaling r and Xj, we shall set 0o — and C3 = 1. 
The remaining parameters can be explicitly given in terms of {Vo,f],a}, 

2(D-2)(q + »y) 2(a + r/) 2 



(a + ??) 2 + 2(D - 2) " ' (a + ??) 2 + 2(D - 2) ' 

°2 = : 7, Q = TTT) 1 — —, > 

a + 77 22 + a z + o;?7 

r K)[(« + r / ) 2 + 2(D-2)] 2 

1 (D - 2) (2 + a 2 + c«7)[2(D - 2) + (D - l)a 2 - r] 2 (D - 3) + 2ar/] ' 1 ' J 

It can be easily seen that there is no scale invariance in such backgrounds. We will 
restrict to 77 > and Vq > without loss of generality, which implies that <fi must diverge 
to positive infinity at the horizon. 

In the specific limit rj = 0, (3 = 2, the scalar potential becomes constant and the scale 
invariance of the solutions can be restored: 

• When we set a — > +00,7 = 2 with vanishing flux, the scalar becomes constant 
and the resulting solution is AdSi>; 
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When we set a — > 0,7 = with flux through R D 2 , <fi becomes constant and the 
resulting solution is AdS 2 x R D ~ 2 ; 

When we set a arbitrary, < 7 < 2 with flux through R D ~ 2 , the scalar (J) ~ log r and 
the resulting solution is the modified Lifshitz solution, whose dynamical exponent 



Before coming to practical calculations we should determine the range of parameters. 
Firstly, by requiring that 0(r) — > +00 for small r and the flux to be real, one can impose 
the bound on a in terms of fixed rj, 

-77 < a < 1]. (2.8) 

V 

Secondly, when the flux is zero, 

a = \~^ /3 = 7= 2 + (D 4 -2)^ " 7 - 

we should require 7' > 1 to ensure a well-defined boundary in the sense of AdS/CFT^, 

2 



if < 



D-2 

Combining constraint derived in the general background in the beginning of this section, 
we can obtain the following complete restrictions 



K/3<2, 7</?<2, -V<a<^~V, < 77 < y ' (2.9) 

We will impose such constraints in the subsequent calculations. 
The near extremal solution can be obtained in a similar way, 

D-2 

_ 4- n n ri 
C x rM(r, 

1 J y ' 8=1 

where 



ds 2 = -Cyf{r)dt 2 + + C 3 r^ £ dxl (2.10) 



/(r) = l-(^±r, u; = j 9- 1 + ^7, (2.11) 
r 2 



1 The reason why such a solution is called "modified Lifshitz" is that the scalar held must be constant 

in Lifshitz background, which is required by scaling symmetry 43!. Properties of the modified Lifshitz 

solutions have been studied in [44] and [20] . 
2 For details see [371 . 
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and the other parameters and fields remain the same as the extremal solutions. One can 
easily get the temperature 

T = -^CiiorJ" 1 , (2-12) 

and the entropy density 

S BH 1 ^ 7 



(2.13) 



3 Thermodynamics of probe D-branes 

In this section we will investigate thermodynamics of massless charge carriers in the 
backgrounds reviewed in previous section. According to AdS/CFT, Nf probe D-branes 
correspond to Nf fields in the fundamental representation of the gauge group in the 
probe limit Nf <C N [HI [46]. An efficient method for evaluating the DC conductivity 
and DC Hall conductivity of probe D-branes was proposed in [17] and [H]. Moreover, a 
holographic model building approach to "strange metallic" phenomenology was initiated 
in [19], where the bulk spacetime was a Lifshitz black hole and the charge carriers were 
described by D-branes. Here we will consider probe D-branes as massless charge carriers 
and explore the thermodynamics in the near-extremal background. 

The dynamics of probe D-branes is described by the Dirac-Born-Infeld (DBI) action 

Sdbi = -N f T D Vol{E) J dtdrd q xe-' t> y/-det{g ab + 2ira'F ab ), 

= -^/«e-V-det( 9 . + 2^), (3.1) 

where To denotes the tension of D-branes, g a b is the induced metric and F a b is the U(l) 
field strength on the worldvolume. In the second line we set r g = ]V/TdVo1(£), where 
Vol(E) denotes volume of the internal space that the D-branes may be wrapping. Further- 
more, we assume that the D-branes are extended along q < D — 2 spatial dimensions of 
the black hole solution. If q < D — 2, the fundamental fields are propagating along certain 
g-dimensional defect. We will introduce a nontrivial worldvolume gauge field A t (r) and 
absorb the factor 27ia' into F a b. Since we are not studying realistic string theories, the 
Wess-Zumino terms will be omitted in the following discussions. 

Before proceeding we should make sure that the backreaction of the probe branes onto 
the background can be neglected. Our discussion is along the line of [49J. Expanding the 
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DBI action to quadric order of F rt in the background, we can obtain 

Sdbi = -T e « J dtdrd q xe^^^l + g u g rr F 2 t , (3.2) 

To avoid backreaction of the probes on the background, the stress energy of the probes 
must be smaller than that generating the bulk spacetime. It can be easily seen that the 
stress energy of the original background~ £p _2 |A|, where £p denotes the Planck length in 
/^-dimensional spacetime and A is the corresponding cosmological constant. Therefore by 
varying the quadric action of the probes with respect to g tt , we can arrive at the following 
condition 



y/1 + g u g rr F r 2 t r cS 
One can see that as long as the effective tension r c fj is sufficiently small, the backreaction 
can be neglected. 

In this background configuration, after performing the trivial integrations on dtd q x and 
dividing out the infinite volume of R 1 ' 9 , we can obtain the action density. 

S=-t cS J drr m \[Y^, (3.4) 

where 

m l r,„ r + + +2(£-2) 

m =2 iq ~ C2= (a + ^ + 2(D-2) ' (3 ' 5) 
and the prime denotes derivative with respect to r. Then the charge density is given by 

_ SC r™A' t 



We can also solve for A' t {r) 



d p 

, — , d = — . 

Vr 2m + d 2 r eS 



K= r^—^ , d = i-. (3.7) 



By plugging in the solution for A' t (r), we can find the on-shell action density 



S'on-shell = — T e ff / dr , ■ (3.8) 

J Vr 2m + d 2 V ; 

Following the methods used in [50J , some interesting physical quantities like the chemical 
potential and the free energy, can be evaluated analytically. We will see that this is still 
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the case for our background. The chemical potential is given by 



poo 

/ A' t dr, 



11. 1 r 2m 



= ^^(s-^ + wt)' (3 ' 9) 

where 

/i = d^So(m), BoH = \b{1 + ^\- 2^)- (3-10) 

Notice that in order to obtain the above results, we have made use of the following useful 
formulae for Beta function and incomplete Beta function 

B(a,b) = = rd«(l + «)-<*»V-\ 

T(a + b) J 

nx/(l-x) 

B(x;a,b) = / du(l + u)- (a+b) u a -\ 

Jo 

as well as for Hypergeometric function 

B(x;a,b) = a~ 1 x a 2F 1 (a, 1 — b; a + 1, x), 

x 

2 F 1 (a,b;c,x) = (1 - x)~ a 2 F 1 (a, c - b; c, — ). 
After choosing the grand-canonical ensemble, the free energy density is given by 

poo r 2m 
^ = — -^cm-shell = r cS / dr =, 

J r+ Vr 2m + d 2 

- ^-ldh3 r2 " +1 > H1 + ^ 2 + ^-P- (3 - U) 

where 

n = - , Teff d 1+ ^B(l + —,--—) = - T f fj,™+\ (3.12) 

2(m + l) v 2m' 2 2m ; (m + l)B (m) m ™ V ; 

Moreover, other thermodynamic quantities can also be calculated from the thermody- 
namic relations. The charge density can be written as 

'=-!r T -<i&> m = T -* (313) 

which is consistent with previous result. The entropy density is given by 

dQ O , 47T x 1 2-/3 T cff . All , 2m+l 2(m + l)-/3 
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Notice that when /3 = 2, there exists a nontrivial contribution to the entropy density at 
T = like those observed in [35] and [36]. On the other hand, the entropy density is 
vanishing at extremality as long as ft ^ 2. The specific heat is 

^ds p(2 — 3) , Air , i 2-0 ri(2m + 2 — 3) , An , 2m+i 2^+2-0 
c v = T— = PK ^ — FiTFi + eff ; ^ (— 7=TT^. 3.15 

It is well known that for a gas of free bosons in q spatial dimensions, the specific heat at 
low temperature is proportional to T p , while for a gas of fermions the low temperature 
specific heat is proportional to T, irrespective of p. When /3 7^ 2 and r c fj is sufficiently 
small, the first term dominates. One can easily obtain /5 = 3/2 when the specific heat is 
proportional to T. Then the parameter a can be expressed in terms of rj 



ai± = (2D - 5)77 ± y/2(D - 2)[2(D - 2)r] 2 - 1], (3.16) 
Combining with (12.9j) . we can arrive at the following conclusions: 

• When 

1 / 2 

< V < 



2(D -2) ~ ' y 3(£>-2) 
both ai+ and aj_ are permitted solutions; 

When 



2 / 2 

< ?7 < 



3(D -2) - ' V D - 2' 
only ai_ is a permitted solution; 

When 

1 



< 77 < 



2(D-2)' 

there is no solution, which means that we cannot realize cy oc T in this regime. 

When f3 = 2,i] = 0, the second term provides the only contribution. The linear depen- 
dence on T fixes m = 1/2. Note that in this limit, 7 = 2a 2 /{a 2 + 2(£> - 2)) = 2/z. So 

-4(D - 2) ± V16(D - 2)2 + gg - 2)(2g - 1) 
" 2± = 2(2^3^ > ( 3 - 17 ) 

By taking into account of (12.91) . we can find that only «2+ is permitted. 
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We can also formally evaluate the "speed of sound". In grand-canonical ensemble, the 
pressure is given by 



p = -O = -. ~^r~~ . , fj£ + \ (3.18) 



while the energy density is 



e = Q + w = ( ™* ( , m ^ +1 - (3.19) 
(m + l)Bo(mj m 

Therefore, 

£ = m p, => c 2 = = _. (3.20) 

cte m 

However, it was emphasized in [3B] that this quantity is only the speed of normal/first 
sound in the relativistic case z = 1. Actually the speed of normal sound is dimensionful 
in a system with z > 1. We will calculate the holographic zero sound in the next section. 



4 The holographic zero sound 

In this section we will calculate the holographic zero sound in the anisotropic background 
at extremality. The basic strategy is to consider fluctuations of the worldvolume gauge 
field on the probe D-branes in the background with nontrivial A t . Such analysis was 
performed for AdS p+ 2 background in [3T] and for Lifshitz background in [36]. We will 
calculate the holographic zero sound in a similar way and classify the behavior of the zero 
sound in different parameter ranges. 

4.1 The retarded Green's functions 

Zero sound should appear as a pole in the density- density retarded two-point function 
Gft((jj, k) at extremality [21] ■ In [SB] the authors provided a general framework for evalu- 
ating the corresponding retarded Green's functions with background metric 

ds 2 = g tt dt 2 + g„dr 2 + g xx ^ dx i- 

i=l 

Here we will take the nontrivial dilaton into account. The symmetries in the spatial 
directions allow us to consider fluctuations of the gauge fields with the following form 

A^r) ->■ A M (r) + a M (t,r,x), 
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where x denotes one of the spatial directions. The quadratic action for the fluctuations 
is given by 

S 2 - f dtdrd q xc-' t> a q/2 \ 9rr ^ x - + \9tt\9rra? , 2) 

where f tx = d t a x —d x a t . Note that we are working in the gauge of a r = 0. After performing 
the Fourier transform 



a,(t,r,x) = J ^e-^a^k), 



the linearized equations of motion can be written as 

— <j> <?/2 1 I / q/2—l 

o r e fl'aia; \9u\9rr a t l e 5Vr /,2 , i \ n //l o\ 

9r[ n7T; " /i i A ^ k at + ^) = °' ( 4 - 3 ) 

U0«l0rr - A ) ' V\9tt\9rr ~ K 

_0 q/2-1, I / _0 9/2-1 

Q r e \9tt \ a x e (7xx g r r / 2 , r \ n //l /l\ 

^nTT^ + /i 1 iTf ^ a * + wA;a ') = °- ( 4 - 4 ) 

U&t I ffrr "A)' Vl^ttlffrr ~ A 

In addition, the following constraint can be obtained by writing a r 's equation of motion 
in a r = gauge 

g rr g xx ua' t + (\g tt \g r r ~ A?)ka' x = 0. (4.5) 

The above equations are not independent, as we can obtain the equation of motion for 
a t by combining the constraint equation and the equation of motion for a x . Therefore 
it is sufficient to solve the constraint and the equation for a t only. By introducing the 
gauge-invariant electric field 

E(r, u, k) = ua x + ka t , 
we can obtain the equation of motion for E 

E" + [d r H e ~?%2 2, )\ E ' - r^(u 2 k 2 - u 2 )E = 0, (4.6) 
u[k l u l -uj l ) \g tt \ 

where 



/ \9tt\9rr I a <-r\ 

U ( r ) = \ — • ( 4 - 7 ) 



9rr9x X 

Moreover, the quadratic action can also be expressed in terms of E, 

/_<£ £=2 1/2 I, 
A A Al 6 9 2 9rr r F 2 . \9tt\ ^2i /, m 

dM u [E + g rr (u^-^) E ] ' (4 ' 8) 
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For our specific background, the metric can be rewritten in terms of the new radial 
coordinate z = 1 jr as follows 

r r/ 3 - 4 i D ~ 2 
ds 2 = -^ d e + Z —dz 2 + — Y dx 2 . (4.9) 

zi 3 Ci z~< ^ 1 v ' 

1 i=l 

In the new coordinate system, the solution of the worldvolume gauge field and the function 
u(z) are given by 

* = z\\ + d 2 z 2m Y U = z^{l + d 2 z 2m Y ^ ' 

where dot denotes derivative with respect to z. Integrating the quadratic action by parts 

* = T fi ^Mf-X ™- <4 - u) 

introducing a cutoff at z = e and taking the limit e — > 0, the quadratic action turns out 
to be 

_ f 2—m 

S a ? = — f- / doudk—EE. (4.12) 

After imposing the incoming boundary condition at the "horizon" z -)• and plugging in 
the solutions of a M , the retarded correlation function reads [51] 

By defining 

Tl(u,k) = J ^S a i, (4.14) 
the retarded correlation functions can be written in terms of II (u , k) 

G? t (io,k) = k 2 Tl(Lo,k), Gg(u,k)=ukU( U ,k), G^ x (u,k)=u 2 U(u;,k). (4.15) 



4.2 Matching the solutions 

In order to evaluate the retarded correlation functions, we should try to solve (14.61) . whose 
analytic solutions are always difficult to find. We will leave the numerical work to section 
5, while here we will obtain the low-frequency behavior of II(u;, k) by solving (14. 6 j) in 
different limits and matching the two solutions in an overlapped regime, following the 
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spirit of [31] and [36] . To be concrete, we will solve ( 14. 6 p in the limit of large z and then 
expand the solution in the small frequency and momentum limit. Next we will take the 
small frequency and momentum limit first and then perform the large z expansion. The 
integration constants can be fixed by matching the two solutions. 

First let us take z — Y <x>, which leads to the following equation for E 

eS-O + Ie+^E-O. (4.16) 
The solution can be given in terms of a Hankel function of the first kind, 



E^D^m), ^c^rr' 1 ' -\-w=ty (4 - 17) 

In the limit of small frequency with v ^ 0, the asymptotic expansion reads 
E = D n T{- + — - — -r^l-ztan 



2 2{f3-l) J y 2(/3 - 1) 



^°r(7T7^-7T " -)( ^ S -^z^-\ (4.18) 



7T v 2(/3-l) 2 jy 2C 1 (r3-l) 

It should be pointed out that the case of v = must be treated separately. In this case 
the corresponding parameter is given by 



—(D - 1)77 + ^/{D - 2) V + 2 - 2). (4.19) 



Now the expansion contains a logarithmic term 

2i 

E~D + -AKMu;/- 1 ) - log(2d(/3 - 1)) + >y E ), (4.20) 

7T 

where 7^ is the Euler constant. 

Next we take uz 13 ^ 2 <C 1, kz 1 ^ 2 <C 1 with cok~^^ being fixed. Then the last term in ( 14.61) 
can be neglected and the equation of E becomes 

E + [d. H e ""fd \ 9ul9 i /2 )]E = 0. (4.21) 
When m/1 and m ^ 7 — /3 + 1, the solution is given by 

7 m - 1 1 1 1113 1 

1 21 1 m K yJ\ + d 2 z 2m m-1 v 2'2 2m 2 2m !> 

m+n 11 q 

-^^—^4* ^ + + ^ -rf 2 ^)] , (4-22) 

m + n 2 2 2m 2 2m J 
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where n = ft — 7 — 1. For either m^lorm^-^ + l, the powers of 2 do not match, 
which will be displayed in Appendix A. We will make use of the following useful formulae 
for the asymptotic expansion 

2 F 1 (a,b;c,x) = (1 - x)~ a ^^ f ~ a ] 2 F l (a, c-b;a-b + l, 1 ^ 



r(6)r(c-a) z 1V ' ' ' 1-x 

_ b r( c )r(a-6) ni/L _ __ L _ 1 1 



+ (1 - x)-" ; yv TT2F 1 (b, c- a; 6 - a + 1 



r(a)r( c -6) z 1V ' ' ' 

2 Fi(-, -; I -x 2 ) = X" 1 log(z + Vl + x 2 ). 
Therefore the large z limit is given as follows when n^O, 

£ ~ L>i + D 2 [^f - ^ - + " - JL)]. (4.23) 

mo rza 2ma 2 2m 2m 

When n = the expansion also contains a logarithmic term 

= °> + B ^ " £ + (^55 ^ - >o g (--)].(4.24) 

To evaluate the retarded correlation functions, we need small z expansion of the solution. 
It can be seen that the second term in the expansion always tends to zero more rapidly, 
being irrespective of n = or not. Therefore we have 

E~ D 1 + D 2 ^-z m ~\ (4.25) 
m — 1 

Assuming that m > lo, the leading order behavior of E reads E(e) — Di,E(e) ~ 
D 2 C\k 2 t m ~ 2 , so the quadratic action turns out to be 

_ r 2—m 

S a 2 = -— dujdk——EE 
2 J k 2 

TeffCl 



dudkDtD 2 , (4.26) 

Z J 

Thus 2 

n(w ' fc) =HjEW Sa2 = jD! Sa2 ^ ( 4 - 27 ) 



'the m < 1 can be dealt with in a similar fashion, see footnote 7 of 136 
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The relation between the integration constants D\ and D2 can be obtained by matching 
the expansions of the solutions in different limits and eliminating the other integration 
constant D . 

Finally we summarize our result for U(u,k) 

When n 7^ the parameters are given by 

r Ci// n 7 1 7 1 

1 " 7^(2^09-1))^ l 2(/3-l) 2 J l 2(/3-l) + 2 J) 

r ggj r / 7 Ixp/ 7 1x^/1 n n_, 

2 " 2 7 rm(2C 1 (/3-l))^ l 2(/3-l) 2 J l 2(/3-l) 2 J l 2 2m' 2m J ' 

5 3 = * + tajl 2 (/-l) ' Go(w) = W ^' (429) 

When n = the parameters are given by 

2CiQ9-l)/*o . . . 1 

<H = , o 2 = 1, o 3 = — , 

irm 7r 

GoM = c 2 log(^ 2 ), 6 = 2CKP-ir (430) 
4.3 Zero sound 

The dispersion relation of the holographic sound mode is given by setting the denominator 
of H(oo, k) to vanish. Similar to the situations discussed in the value of n determines 
which term in the denominator dominates. 

It can be seen that when n < 0, the u 2 term dominates. Then we can expand k(uj) 

k = ±uM[l + #uT 1+ ^ + 0(uT 2+ ^)], (4.31) 
V 0i 2d 2 



and then invert to find 



u(k) = ±k x ^ - -f (-±)*&=Vkv=i + <D(k- i+ —i). (4.32) 

Notice that when n < 0, p — 1 < 7, so at low momentum k?- 1 < k, that is, the real part is 
bigger than the imaginary part. Therefore this mode describes a quasi-particle excitation. 
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As a check of consistency, we can take the specific limit = 7 = 2, C\ = 1, A = 0, m = q, 
which just gives the AdSn background. One can obtain 

which agrees with [3T]. The speed of the holographic zero sound is given by 



v l= 5 -l = 9l d ^ r (lk) r (l h) {A m 

d2 m 1 (2 + 2^) L y~^> 

It can be seen that in the relativistic case = 7 = 2, the speed of zero sound coincides 
with the speed of normal/first sound. One special case is m + n = 1, where all the T 
functions cancel and the speed of zero sound reads 

v 2 = 9l d ^^ (4.35) 
m 

which turns out to be finite as n — > 0. When m + n > 1 as well as n < 0, r(— ^-) has a 
pole as n — > 0, so f q goes to zero from above. 

When n > 0, Go( w ) dominates, then 

/ O'i. \ — 1 tt- (Jo n 2n 

fc(w) = ±(^) 2 u 1 '^^) (1 + -^c^ + C(o;^)). (4.36) 
Oi 2d 3 

Inverting the relation above, 

, 8^ s P- 1 g(g-l) 8 — 1 5-} / 5l s 2(/3-l)-n 2(g-l) + 4n 

w = ay»-i;-njfe 3 (/j-i)-» + ^: — ^(J.)2W-D-"fc2 W -i)-™ + a</>-i)-») (4.37) 

ids 2(0-l)-nd 3 K d 3 / 

Notice that since 8\ is real and 5 3 is complex, the leading term has a complex coefficient. 
Furthermore, the real and imaginary parts are of the same order, hence this mode is not 
a quasi-particle. 

Finally when n — 0, 



k = ±^=^5 2 + \og(5u 2 ). (4.38) 



Expanding for small u, 



k(cu) = ±^=y/6 3 ]og(6u*) - ^(5 3 \og(5cu 2 ))^ + 0{uj\og-V\5uj 2 )). (4.39) 



It can be seen that the dispersion relation differs from the holographic zero sound mode 
by logarithmic terms. 
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5 A numerical survey of Fermi surface 



In previous section we observed a sound-like excitation in the regime of n < 0. It is known 
that such zero sound mode is associated with the deformation of the Fermi surface away 
from the spherical shape. The theory of normal Fermi liquids tells us that the jump in 
the distribution function can be observed as a singularity in the retarded current-current 
Green's function in the u = limit. To see if we can observe such Fermi surface we need 
the complete solution to the equation for the gauge fluctuations (I4.6p with co = 0, at least 
numerically. 

We investigated thermodynamics of probe D-branes in section 3, where we found that 
when (3 2, the specific heat was proportional to the temperature T under certain 
conditions, which is just the behavior of Fermi liquids. As we have many groups of pa- 
rameters (a, 77) which lead to specific heat linear in T, we choose the following parameters 
a = 3/2,77 = 1/2 in .0 = 4 dimensional spacetime as an example. For simplicity we also 
fix q = 1 and Vb = 1. Then the equation for the gauge fluctuations ( 14. 6 j) in the u = 
limit is reduced to 

„ 3(r 3 -2d 2 ) , I5k 2 y/r 
E + 2r{r* + <P) E ~ 8(r* + <P) E ~ °' (5 - 1} 
In the near horizon region r — > 0, the perturbative solution of the above equation is given 
by 

E = ci + c 2 r 4 . (5.2) 

Here, we choose the boundary condition E = at the horizon, which is compatible to 
impose the incoming boundary condition for w^O. So we set c\ = 0. At the boundary 
r = 00, the asymptotic solution of (15.11) becomes 

E = A + Br- 1/2 . (5.3) 

Notice that since the first term is finite at the boundary, it corresponds to the source term 
and the coefficient of the second term implies the vev of the dual gauge operator. Then, 
the Green function is proportional to —B/A. In Figure 1, we plot the dependence of the 
Green's function on the momentum k for d — 1. 

From the figure it can be easily seen that the characteristic structure does not appear 
in a wide range of k, while the specific heat and zero sound exhibit typical features of 
Fermi liquids. Such phenomenon was also observed in [33J . It was pointed out in [33j that 
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Figure 5.1: Ratio B/A as a function of k. 

one difficulty in applying Landau's theory of Fermi liquids was the assumption that the 
particle number should be conserved as the strength of the interaction was varied, which 
was not obvious for the case they discussed. Their conclusion seems to be still applicable 
to the present case. 



6 AC conductivity 

In this section we calculate the AC conductivity by making use of the correlation functions, 
which can be seen as a by-product of section 4. Furthermore, we will take the limit of 
zero density and compare the results with those appeared in previous literatures. 

It can be easily seen that the current-current correlation function is given by 

k = 0) = c 2 n(c, k = 0) cc - 7^r] y (6- 1 ) 
Therefore we can obtain the following expressions in the small frequency limit, 

G* x (u ->■ 0, k = 0) oc -TeffC^w^ n > 0, (6.2) 



G*(u -> 0, k = 0) oc - T f°lf n = 0, (6.3) 
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G* x (u -> 0, k = 0) oc -r eff Ci52 \ n<0 (6.4) 
Recalling that the definition of AC conductivity is given by 

a(uj) = --Gt(u^0,k = 0), (6.5) 
ui 

we can arrive at the following results 

1 — 1 

<t(o;) oc iT e ffCi<5^ a;^ 31 " , n > 0, (6.6) 

<7(w) OC ITeffCitfg ^^(log^W 2 ))- 1 , 71 = 0, (6.7) 

<t(u) oc zreffCi^ o; -1 , n < 0. (6.8) 

As a check of consistency, we consider the specific limit, = 2,7 = 2/z,n = 1 — 2/z. When 
n > 0, z > 2, cr(u;) oc ZT e fjC 1 (^ 1 u;~ 2 ' z , which agrees with the result obtained in [36]. The 
behavior of the AC conductivity is qualitatively similar to that investigated in [36] . To be 
concrete, when n 7^ 62 is real while £3 is complex. The conductivity is purely imaginary 
when n < and has a simple pole at zero frequency. Then according to Kramers-Kronig 
relation, one can conclude that the real part of the conductivity, and hence the spectral 
function, consists of a delta function at zero frequency. When n > the conductivity, and 
hence the spectral function has a power-law dependence, and no explicit quasi-particle 
excitation exists. 



On the other hand, it was observed in [52] that after transforming the equation of a x 
into a Schrodinger-like form, the AC conductivity was directly related to the reflection 
amplitude for scattering off the potential. For our system, it can be observed that the 
equation for a x at k = becomes 

3/2-1 1 I / 3/2-1 

drl n + 1 ^ = = (6.9) 

r[ (\g tt \9rr-A?y/2 l y/\g tt \grr-A? 



in the original coordinate 



2 D-2 



i=l 
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The asymptotic behavior of a x is given by 

a» = + . J f ] x1 . r^, (6.10) 

r cff (7 - /3 - m + 1) 

where we have introduced a background electric field E x {t) = ReE x (u)e~ luJt , The above 
equation can be put in a Schrodinger-like form, 

d 2 V 



ds 2 



+ u(s)y = n 2 ^, (6.11; 



where 



-7/2 A A ,. 

A = Vr 2m + d 2 , -f=r- /3 -f, n=^-, (6.12) 



and 



* A 1 / 2 v ' ' dr ds' C7i : 

tf( a ) = r 7/2 A -V2 [ l^. (r -7/2 A -l/2^ ) _ ^.(7-7/2-^1/2^^ ( g 13) 

Generally speaking, the AC conductivity can be obtained by numerical methods. However, 
in the zero density limit 

d = 0, X = r m , (6.14) 

one can find 

TT( a \ Uo TT (m- 7 )(m-7 + 2/3-2) 

U(s) = ^, U = , (6.15) 

so the potential possesses the same form as those investigated in [17] and [20]. Assuming 
that the solution asymptotes to AdSo, we can obtain 

Re(a)~w K , « = 2i/b-l, u 2 = U + - (6.16) 

following their approach. In particular, when q = 2, C2 = 0, the potential vanishes and 
the conductivity is constant at all temperatures 

a(u) = t cS = (T , (6.17) 

which agrees with the analysis performed in [53J. 



7 Summary and discussion 

In this paper we explore the zero sound in D-dimensional effective holographic theories, 
whose bulk fields include the graviton g^ u , the U(l) gauge field and the scalar field 
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(f). The solutions possess anisotropic scaling symmetry and they reduce to previously 
known examples, such as AdSp, AdS2 x R^ -2 and "modified" Lifshitz solutions under 
certain conditions. We consider thermodynamics of massless probe D-branes in the near- 
extremal background and clarify the conditions under which the specific heat is linear 
in the temperature, which is a characteristic feature of Fermi liquids. Subsequently we 
study the zero sound mode by considering the fluctuations of the worldvolume gauge 
fields on the probe D-branes. Rather than analytically solving the equations of motion, we 
obtain the low-frequency behavior by solving the equation in two different limits and then 
matching the two solutions in a regime where the limits overlap, following (31] and [36] . 
The resulting behavior of the zero sound looks similar to that investigated in [36], that 
is, when the parameter n = — 7 — 1 < 0, the dispersion relation reveals a quasi- 
particle excitation; while the zero sound is not a well-defined quasi-particle when n > 0. 
Furthermore, we plot the correlation function in D = 4 at u = with fixed parameters 
which lead to linear specific heat. The result is that we cannot observe any characteristic 
structure of Fermi liquids in a wide range of k, which is similar to what was found in [33] . 
As a by-product, we also evaluate the AC conductivity via the current-current correlation 
function, which reduces to previously known results at specific limits. 

By now there are mainly two approaches for studying condensed matter physics in the 
context of holography. One approach can be thought of as "top-down", that is, we 
consider certain exact solutions or brane configurations in string/M theory which possess 
the desired properties of condensed matter systems. The main advantage is that we 
have clear understanding about the dual field theories, while it is difficult to find such 
solutions in string/M theory. A complementary approach can be seen as "bottom-up", 
which means that we consider certain toy models of gravity which possess solutions with 
the desired properties. It allows a parametrization of large classes of IR dynamics and 
provides useful information in the dual field theory side. However, the main disadvantage 
is that the embeddings of such toy models into string/M theory are not obvious, thus many 
things in the field theory side remain unknown. However, the "bottom-up" approach is 
an efficient tool for investigating the AdS/CMT correspondence. 

The background solutions we studied in this paper are exact solutions of theories with 
domain wall vacua. In [37J the author also constructed interpolates between two exact 
solutions of the single- exponential, domain wall gravity theory, which lead to the argu- 
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ment that the domain wall/QFT correspondence [54] can be considered as an effective 
holographic tool which is applicable in settings beyond the regime of domain wall su- 
pergravities. Therefore the domain wall/QFT correspondence can also be taken as one 
specific class of effective holographic theories. Moreover, the author of [37] argued that 
even when the UV completion of some bulk theory was unknown, if the theory admitted 
an approximate domain wall solution at some intermediate value of r then one could use 
domain wall/QFT correspondence to develop a holographic map. Thus it is interesting to 
develop the AdS/CMT holography by making use of this domain wall/QFT correspon- 
dence. In particular, we can establish precision holography in this anisotropic background 
along the line of [55] and study the fermionic correlation functions following [56| [57J |58] . 
We leave such fascinating projects in the future. 
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A Asymptotic expansions of E(z) at specific values 

In this appendix we show that when the parameter m takes some specific values, the large 
z expansions of the solutions E(z) to equation (I4.2ip cannot match that of (I4.17p . Firstly, 
when m = 1 but m ^ —n, the solution to (I4.2ip is given by 

E{z) = D 1 + D 2 [C 1 k 2 ( 1 j9 +logz-log2(l + Vl + d 2 z 2 )) 

V 1 + d A z A 

9 I Z 13 ' 1 „ , 1/3—7 ft — 7 ,9 9n 

2 + -<<*))]■ (A.i) 

2 + /3-7 v 2' 2 ' 2 ' !n K ' 

When performing the expansion, we just focus on the powers of z. The Hypergeometric 
function gives 

~ const + z p ~ 1+l + z~ 2 \ (A.2) 
which does not match that of ( 14.171) . 
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Secondly, when m/1 but m = — 



n ^ 0, the solution is 



E(z) 



D 1 + D 2 [C 1 k 




m ~ l , 1 1 „,1 1 13 1 




1 



+ (1 - (3 + 7) log z - log(2 + 2y/l + cPz^ 1 ^))}. (A.3) 



1-/3 + 7 



The u 2 term gives 



which does not match that of (I4.2ip either. Therefore we just considered the cases with 
m/1 and m 7^ 7 — + 1 in the main text. 
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